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The results of experimental studies in selecting optimum characteristics for diamond drills and optimum treat¬ 
ment procedures for drilling holes 1-70 mm in diameter in hard brittle nonmetallic materials such as quartz, 
ceramics, glass, and glass ceramics are considered. To obtain experimental-mathematical models, a random¬ 
ized experimental design was used at the screening stage, i.e., a half-replica of the full factorial experiment 2 3 , 
and an orthogonal compositional design of the second order of type 3 s was used at the main stage. The dia¬ 
mond drilling process was optimized using simplex methods. The correlation of the estimated values with ex¬ 
perimental data showed the adequacy of the results obtained. 


The use of solid and brittle nonmetallic materials (quartz, 
glass ceramics, glass, and ceramics) calls for upgrade of me¬ 
chanical treatment methods. For this purpose tools with high 
abrasive characteristics based on natural or synthesized dia¬ 
monds are developed and actively applied in the manufactur¬ 
ing industry. The efficiency and quality of diamond treat¬ 
ment of solid and brittle nonmetallic materials can be im¬ 
proved by optimization of treatment regimes. An important 
instrument for solving this problem is a mathematical experi¬ 
mental design using general principles of probability theory 
and mathematical statistics. 

It is advisable to apply probabilistic-statistical methods 
to study complex technological processes and experimental 
data processing to determine the reliability, accuracy, and 
sufficiency of data for making certain decisions and also to 
develop mathematical-statistical models to be used for the 
process optimization. 

The process considered is a complex system, and we will 
attempt to account for the mutual influence of factors to ob¬ 
tain an integrated representation of treatment in a mathemati¬ 
cal form. 

The factors that have an effect on the output parameters 
of the model were identified based on empirical data and re¬ 
sults of analytical modeling [1]. 

The model parameters are the following: 

- axial cutting force P f ; 

- roughness of the surface treated R a ; 

- average width of chips at the hole inlet 
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- specific consumption of diamonds q\ 

- resistance of the drill l. 

The factors influence the model parameters are as follows: 
Xj) cutting velocity v ; 
x 2 ) feed S ; 

x 3 ) size of diamond grains G; 

x 4 ) concentration of diamonds C; 

x 5 ) sort of diamond; 

x 6 ) type of binder material; 

x 7 ) material of the part treated; 

x H ) type of lubricant-coolant (LC) T [ C ; 

x 9 ) lubricant-coolant pressure P [ C ; 

x 10 ) lubricant-coolant flow rate Q LC . 

The total number of factors is too high for varying, and 
consequently studies were carried out to simplify the model 
and reduce the number of its factors. 

A first series of experiment was carried out using a ran¬ 
domized experimental design for a preliminary investigation 
of the response surface on a chosen site. To obtain an experi¬ 
ment design, the factors investigated were split into two 
groups, and a semireplica of a full factorial experiment of 
type 2 5 was used for each group. The initial data are listed in 
Table 1 and the results are listed in Table 2. 

Based on the data from Table 2, scattering diagrams were 
constructed for each of 10 factors, and factors significant for 
the diamond drilling process were identified using the r|-rule 
[2], An example of a scattering diagram is shown in Fig. 1. 
The significant factors of the process are unequally distri¬ 
buted among the controlling variables. As a consequence of 
ranking, the following data were obtained. 
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TABLE 1 


No. 

jCj, m/sec 

o 2 , mm/min 

x v fim 

x 4 , % 

x s 

*6 

x i 

x 8 

x 9 , MPa 

x 10 , liter/min 

1 

3 

10 

400/315 

150 

AS50 

Ml 

22KhS 

EMUS 

0.3 

8 

2 

1 

30 

400/315 

50 

AS 160 

Ml 

K8 

EMUS 

0.1 

8 

3 

1 

30 

400/315 

150 

AS50 

Ml 

K8 

EMUS 

0.3 

2 

4 

1 

10 

400/315 

150 

AS 160 

Ml 

K8 

NGL-205 

0.3 

8 

5 

3 

30 

400/315 

150 

AS 160 

Ml 

22KhS 

NGL-205 

0.3 

2 

6 

1 

10 

200/160 

150 

AS50 

OPM 

22KhS 

EMUS 

0.3 

2 

7 

1 

30 

200/160 

50 

AS50 

OPM 

K8 

NGL-205 

0.3 

2 

8 

3 

30 

200/160 

150 

AS50 

OPM 

22KhS 

NGL-205 

0.3 

8 

9 

1 

10 

200/160 

50 

AS 160 

Ml 

22KhS 

EMUS 

0.1 

2 

10 

1 

10 

400/315 

50 

AS50 

OPM 

22KhS 

EMUS 

0.1 

8 

11 

3 

10 

200/160 

50 

AS50 

OPM 

22KhS 

NGL-205 

0.1 

2 

12 

3 

30 

400/315 

50 

AS50 

OPM 

K8 

EMUS 

0.3 

8 

13 

3 

10 

400/315 

50 

AS 160 

Ml 

K8 

NGL-205 

0.1 

2 

14 

3 

10 

200/160 

150 

AS 160 

OPM 

K8 

NGL-205 

0.1 

8 

15 

1 

30 

200/160 

150 

AS 160 

Ml 

22KhS 

NGL-205 

0.1 

8 

16 

3 

30 

200/160 

50 

AS 160 

OPM 

K8 

EMUS 

0.1 

2 


TABLE 2 


No . 

P f , N 

R a , |im 

mm 

q , mg/cm 

/, mm 

1 

163 

0.9 

0.12 

0.20 

25 

2 

144 

2.4 

0.21 

0.15 

500 

3 

126 

2.1 

0.20 

0.14 

450 

4 

60 

1.5 

0.18 

0.13 

700 

5 

275 

0.9 

0.11 

0.18 

15 

6 

213 

0.6 

0.08 

0.22 

30 

7 

162 

1.9 

0.17 

0.13 

350 

8 

312 

0.8 

0.11 

0.25 

5 

9 

187 

0.6 

0.07 

0.19 

15 

10 

212 

1.4 

0.09 

0.21 

10 

11 

200 

0.8 

0.07 

0.19 

15 

12 

138 

2.4 

0.20 

0.15 

550 

13 

84 

2.3 

0.24 

0.16 

600 

14 

108 

1.5 

0.15 

0.14 

570 

15 

300 

0.7 

0.12 

0.20 

10 

16 

168 

1.7 

0.19 

0.12 

440 


Significant Factors of the Process 


Model parameters 

Axial cutting force . . . 


Roughness of surface treated . 


Average width of chips 
at the hole inlet. . . 


Specific consumption 
of diamond .... 


Drill resistance 


Factors influencing the parameter 

Feed, diamond grain size, diamond 
concentration, sort of diamond, 
material of part treated 

Cutting velocity, feed, diamond 
grain size, diamond concentra¬ 
tion, material of part treated 

Feed, diamond grain size, diamond 
concentration, material of part 
treated, LC flow rate 

Cutting velocity, feed, sort of dia¬ 
mond, type of binder, material 
of part treated 

Cutting velocity, feed, sort of dia¬ 
mond, material of part treated, 
LC flow rate 


P f , N 



Fig. 1. Scattering diagram of P f . 


Based on analytical modeling [1], the following quantita¬ 
tive parameters were chosen to characterize, respectively, the 
material of a product treated, the binder type, and the sort of 
diamond: the microhardness of product material //; the bend¬ 
ing strength of binder a ^ m , and the compression strength of 
diamond powder a . 

Thus, it is necessary to determine the following 
dependences to implement optimization: 

P f =f(S,Qo?\H); 

R a =/(v, 5, G, C, H); 


^ = f(S, G, C, H, Q lc ); 


dia bin 


q=f(v,S,C™,a 


H); 


l=f(y, Q lc ). 
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TABLE 3 


Factor 

Variation 

interval 


Level 


main 

upper 

lower 

Xj, m/sec 

i 

2 

3 

i 

x 2 , mm/min 

10 

20 

30 

10 

x 3 , jam 

100 

180 

280 

80 



(200/160) 

(315/250) 

(100/80) 

x 4 , % 

50 

100 

150 

50 

+5 

30 

56 (AS60) 

86 (AS 100) 

26 (AS20) 

+6 

75 

130 (OPM) 

205 (M) 

55 (Ml) 

x 7 

500 

1000 

1500 

500 



(STM-1) 

(TsM332) 

(sheet glass) 

x 9 , MPa 

0.1 

0.2 

0.3 

0.1 

x 10 , liter/min 

3 

5 

8 

2 


It is known a priori that the dependences considered can 
be described by a polynomial of the second order: 



Fig. 2. Results of modeling the mutual effect of parameters that in¬ 
fluence specific consumption of diamond q depending on the cutting 
velocity v. 


Y = B 0 + 5, Xj + B 2 x 2 + B 3 x 3 + B 4 x 4 + B 5 x 5 + 
B xl x x x 2 + 5 13 X[ x 3 + B xa X\ x 4 + B l5 x l x 5 + B 23 x 2 x 3 + 
B 2 4 x 2 x 4 + B 25 x 2 x 5 + B 34 x 3 x 4 + B 35 x 3 x 5 + B 45 x 4 x 5 + 
B n X\ + B 22 x 2 + B 33 x 3 + B 44 x x +S 55 X!, 


adequacy of the model obtained using the /•'-criterion (Fisher 
criterion), mathematical models of the processes investigated 
were constructed. For instance, such a model for P { has the 
following form: 


where /f are the model coefficients; x ( are the process factors. 

An orthogonal compositional second-order design of 
type 3 5 was accepted as the experimental design. The factor 
variation intervals were determined as a consequence of pre¬ 
liminary experiments and analytical modeling (Table 3). 

In view of the orthogonality of the design, all regression 
coefficients are calculated independently of each other using 
the following formulas: 

i ”o 

h o= IX; 

n 0 u = 1 
N 

k = j =i 


j =i 
N 

j =i 


7 = 1 

where n 0 is the number of experiments at the center point of 
the design; u is the consecutive number of the replicate ex¬ 
periment; Y 0u and Y. are the values of the parameter consi¬ 
dered in the «th and the yth experiments; j is the consecutive 
number of the experiment; i is the consecutive number of the 
factor; x.. is the ith factor in the yth experiment. 

Based on the modeling results, regression equations were 
obtained for the parameters considered. After verifying the 


P { = 22.6710-9.271 x 10 - 3 S+9.69 x 

3.1826 x 10“ 2 G+ 1.62 x 10^ 5 G 2 -2.8056 x 10“ 2 C + 
1.77 x 10" 4 C 2 - 0.42256af a + 2.97 x 10- 3 of a2 - 

3.992 x 10 -3 77+ 7.856 x 10- 6 77 2 + 6.8x 10' 5 5G + 

5.98 x 10 5 SC + 2.73 x lO^^of 2 + 5.485 x 10 ~ 4 SH- 
2.01 x 10“ 5 G5” + 2.45 x 10“ 5 Gof a + 9.24 x 10 - 1 GH + 
2.35 x 10“ 5 Cof a + 3.28 x 10“ 7 CH-2.35 x 10- 6 /7o 4ia . 

Mathematical models for other parameters were calcu¬ 
lated similarly. 

The second stage of optimization of the diamond drilling 
process is identifying the coordinates of the local optima of 
the models constructed using the simplex method. The effi¬ 
ciency of drilling was taken as the criterion of optimality. 
The principal constructions were implemented in the v-S 
coordinate system, and the dependence obtained were used 
to find the limiting surfaces. Plots (surfaces) describing the 
mutual effect of the parameters were constructed for the 
dependences obtained. An example of a set of surfaces for 
specific consumption of diamonds is shown in Fig. 2, which 
represents the limitation surface for a maximum specific 
consumption and the surfaces for the calculation of optimum 
efficiency parameters [1 ) based on v - a * a ; 2) based on 
v-H-3) based onv-cj 1 ”]. 

The analytical and experimental mathematical models 
were used to develop algorithms and specialized software to 
calculate tool parameters and treatment conditions taking 
into account the following limitations: R a < 1.6 pm, 
^ < 0.10 mm, q < 0.30 mg/mm, l > 200 mm. 
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TABLE 4 


Parameter* 

Horseshoe-shaped drill 
of diameter 1-5 mm drilling 

Circular drills 

of diameter 6-70 mm drilling 

glass ceramics 

quartz 

glass 

ceramics 

Cutting velocity, m/sec 

5 

5 

3 

3 

Feed, mm/min 

5 

5 

29 

8 

Size of diamond grains, pm 

100/80 

125/100 

160/125 

200/160 

Concentration of diamonds, % 

350 

350 

250 

250 

Sort of diamond 

AS50 

AS65 

AS50 

AS80 

Lubricant-coolant flow rate, liter/min 

8 

8 

6 

6 


* Binder M was used in all cases. 


Computerized calculations using the models described 
produced optimum values for the tool parameters and treat¬ 
ment conditions (Table 4). 

The axial cutting force calculated for the specified drill 
parameters and treatment regimes is equal to 53 and 67 N for 
drilling holes 3 mm in diameter in glass ceramics and quartz, 
respectively, and 141 and 253 N for drilling holes 26 mm in 
diameter in glass and ceramics, respectively. 

A comparison of the calculated values with experimental 
data shows the adequacy of the results obtained. Thus, the 
axial cutting force (which is the main parameter used in 
adaptive control systems for diamond drilling) in verification 
experiments amounts to 50 and 70 N in drilling 3-mm hole in 
glass ceramic STM-1 and in quartz, respectively, and 130 
and 250 N, respectively, in drilling 26-mm hole in sheet glass 


and in ceramic TsM-332. The mean error in the above expe¬ 
riments was 5.5%. 

The use of the modeling results in the form of specia¬ 
lized software for the calculation of efficient cutting regimes 
for various combinations of tools and materials makes it pos¬ 
sible to improve the process efficiency, to reduce the number 
of tool dressings, and to decreases the probability of emer¬ 
gency failures. 
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